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Scientific  Progress 

Our  specific  aims  involved  the  development  of  accurate  sharp  interface  methods  based  on  Cartesian  grids  and  immersed 
interface  treatment  for  incompressible  flows  with  a  free  boundary  or  a  moving  interface. 

One  of  the  most  remarkable  progresses  is  that  we  have  developed  an  efficient  augmented  IIM  for  modeling  the  fluid  structure 
interactions  between  a  fluid  flow  modeled  by  the  Stokes  equations  and  a  porous  media  modeled  by  the  Darcy’s  law.  There  are 
many  applications  of  the  coupled  Stokes  and  Darcy  flow,  for  example,  the  underground  water  flow,  oil  recovery,  fluid  flow 
through  tissues  and  cells.  Our  new  idea  is  to  introduce  several  augmented  variables  along  the  interface  between  the  fluid  flow 
and  the  porous  media  so  that  the  problem  can  be  decoupled  as  several  Poisson  equations.  The  augmented  variables  are 
chosen  so  that  the  Beavers-Joseph-Saffman  (BJS)  and  other  interface  conditions  are  satisfied.  In  the  discretization,  the 
augmented  variables  have  co-dimension  one  compared  with  the  primitive  variables  and  are  solved  through  the  Schur 
complement  system.  Our  work  is  the  first  using  Cartesian  meshes  and  fast  Poisson  solvers  fur  such  a  coupling.  The  method 
can  easily  apply  to  other  coupled  systems  and  applications. 

Working  with  Dr.  Jianlin  Xia  of  Purdue  University,  we  have  developed  an  efficient  preconditioning  technique  for  the  Schur 
complement  system  arises  from  the  augmented  IIM.  The  preconditioning  techniques  have  shown  to  be  very  efficient  for  all  the 
application  problems  that  we  have  tested  including  the  flow  past  obstacles,  contact  line  problems;  drop  spreading,  and  the 
couple  of  the  Stokes  flow  and  Darcy’s  law.  Equipped  with  the  new  preconditioning  technique,  the  Augmented  IIM  has  become 
a  direct  challenger  of  the  traditional  boundary  integral  method. 

Furthermore,  we  continue  our  research  in  the  immersed  finite  element  method  (IFEM).  We  have  developed  the  Robin-Robin 
domain  decomposition  method  for  Stokes-Darcy  coupling  using  modified  Cartesian  meshes.  Based  on  Z.  Li’s  early  non- 
conforming  IFE  method,  we  have  developed  a  symmetric  and  consistent  IFE  method  which  is  second  order  accurate  for  elliptic 
interface  problems  with  homogeneous  and  non-homogeneous  jump  conditions.  This  is  a  significant  step  forward  in  the  IFEM 
research.  The  paper  has  been  accepted  by  the  Journal  of  Scientific  Computing.  Our  rigorous  proof  of  the  immersed  boundary  is 
first  order  accurate  in  the  infinity  norm  for  elliptic  interface  problems  has  been  accepted  by  the  journal  of  Mathematics  of 
Computations.  This  may  be  a  40-year  unsolved  problem. 

We  have  developed  the  direct  sampling  method  for  inverse  medium  problems, 
specifically  determine  distribution  of  scatters  by  near  field  scattering  data 

from  a  single  plane  incident.  We  developed  the  probing  index  function  over  the  medium  domain  that  provides  an  estimate  of  the 
distribution  of  scatters. 

We  use  the  fundamental  solution  to  Helmholtz  equations,  i.e.,  the  probing 
index  uses  the  inner  product  of  scattering  data  with  the  fundamental  solution 

at  the  measurement  boundary.  The  effectiveness  of  the  probing  function  is  tested  for  several  simulated  data.  We  also  extend 
the  probing  method  for 

the  EM  inverse  scattering  problems  based  on  Maxwell  equations. 

We  have  developed  a  comprehensive  Lagrange  multiplier  theory  for 

a  general  class  of  nonsmooth  and  non-convex  optimizations  including  Non-Newtonian  mechanics  and  sparsity  optimizations. 
Nonsmooth  solution  method 

is  increasingly  important  to  construct  the  sharp  and  enhanced  solutions  for  scientific  computations  in  general.  Based  on  the 
theory  we  developed  a  semi-smooth  Newton  method  in  the  form  of  primal  dual  active  set  method.  We  have  demonstrated  the 
effectiveness  and  applicability  of  our  proposed  algorithm  for 

tomography  problems  and  control  problems  and  nonlinear  least  square  solution  for  reconstruction  of  medium  coefficients. 


We  have  developed  the  weak  formulation  for  the  fluid-structure  interaction 

based  on  the  incompressible  Navier  Stokes  equations.  The  formulation  provides  a  systematic  and  accurate  model  for  the 
iteration.  Based  on  the  formulation 

we  introduced  an  energy  preserving  weak  solution  and  we  developed  the  numerical  simulation  method  for  a  general  class  of  the 
fluid-structure  interaction  problems.  We  have  a  very  promising  preliminary  result  and  more  numerical  tests  are  currently  carried 
out. 
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